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Abstract 

Precipitation kinetics of AbjZr and AI3SC in aluminum supersaturated solid solutions 
is studied using cluster dynamics, a mesoscopic modeling technique which describes 
the various stages of homogeneous precipitation by a single set of rate equations. 
The only parameters needed are the interface free energy and the diffusion coefficient 
which are deduced from an atomic model previously developed to study the same 
alloys. A comparison with kinetic Monte Carlo simulations based on the vacancy 
diffusion mechanism shows that cluster dynamics correctly predicts the precipitation 
kinetics provided a size dependent interface free energy is used. It also manages to 
reproduce reasonably well existing experimental data. 
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1 Introduction 



Transition elements are added to aluminum alloys so as to obtain small ordered 
precipitates and by this way decrease the alloy sensitivity to recrystallization. 
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Among these elements, zirconium and scandium are the most efficient ones. 
They lead to the precipitation of respectively Al 3 Zr and Al 3 Sc with the Ll 2 
structure [1-6]. As recrystallization is controlled by the precipitate density 
and size, it is important to model precipitation kinetics in Al-Zr and Al-Sc 
alloys, in order to optimize the heat treatment and better control the alloy 
final state. 

In a previous study [7,8], one of us built an atomic model for the above two 
binary alloys and studied precipitation kinetics using Monte Carlo simulations 
based on the vacancy diffusion mechanism. The drawback of such an approach 
is that only short annealing times can be simulated which implies using high 
supersaturations. We showed that the classical nucleation theory manages to 
reproduce atomic simulations giving thus good confidence in the ability of 
mesoscopic models to quantitatively predict the nucleation stage of precipita- 
tion kinetics. In the present work, starting from the same atomic model, we 
build a mesoscopic modeling, the cluster dynamics [9-11] which is based on 
a set of rate equations describing in the same framework the three different 
stages of precipitation, i.e. nucleation, growth and coarsening. It requires only 
a limited number of parameters, the interface free energy and the solute diffu- 
sion coefficients, which can be quite easily deduced from atomic parameters. 
The range of supersaturations and annealing times that can be simulated are 
thus extended allowing a comparison with experimental data. Cluster dynam- 
ics technique has already been used to model precipitation kinetics in Al-Zr 
and Al-Sc alloys in the spirit of nucleation theory [12]. The way we use clus- 
ter dynamics in the present paper is distinct, as will be shown below, and is 
based on the description of the alloy as a gas of clusters. Unlike the nucleation 
theory, our approach leaves no arbitrary choice for the parameters used. 

In the first part of this article, cluster dynamics modeling is presented. We 
show then how precipitation kinetics obtained with this method compare with 
kinetic Monte Carlo simulations for both Al-Zr and Al-Sc systems and try to 
improve the kinetic description so as to obtain a better agreement. Finally, 
available experimental data are compared with simulations. 



2 Cluster dynamics 

In its strict sense, cluster dynamics rests on the description of the alloy under- 
going phase separation as a gas of solute clusters which exchange solute atoms 
by single atom diffusion [9-11]. Clusters are assumed to be spherical and are 
described by a single parameter, their size or the number nx of solute atoms 
X they contain. In such a description, there is no precisely defined distinc- 
tion between the solid solution on the one hand and the precipitates on the 
other hand, at variance with the classical nucleation theory: the distribution 
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of cluster sizes is the only quantity of interest (for a detailed discussion see 
Ref. 13). In an undersaturated solid solution of nominal concentration z| a t 
equilibrium, the cluster size distribution is given by (cf. appendix A) 

C nx (4) = exp [- (G nx - 2n x /i(4)) I^t] , (1) 

G nx being the free energy of a cluster containing nx atoms and /J-(x x ) = 
(Hx(xx) ~ /^aiC^x)) /2 the effective chemical potential in the solid solution, 
i.e. half the difference between the chemical potentials of the solute and the 
solvent atoms. In a supersaturated system, the cluster size distribution evolves 
in time as discussed below. 



2.1 Master equation 



The cluster dynamics technique describes the precipitation kinetics thanks to 
a master equation giving the time evolution of the cluster size distribution 
[9-11]. When only monomers can migrate, which is the case for Al-Zr as well 
as Al-Sc systems [7], the probability C nx to observe a cluster containing n x 
solute atoms obeys the differential equations 

^n x -l-m x — ■Ai x ->n x +l> ^ n X > 2 (2a) 

— 2Ji_ >2 — ^2 Jn x ->n x +l, (2b) 
n x >2 

where J nx ^nx+i is the cluster flux from the class of size nx to the class nx + 1. 
This flux can be written 

J-rix^nx+l Pnx^nx ^nx+lCrix+1 ; (3) 

Pn x being the probability per unit time for one solute atom to impinge on a 
cluster of size nx and a nx for one atom to leave a cluster of size %• 

In order to handle large cluster sizes, the set of differential equations 2 and 3 
are integrated by considering cluster sizes in a discrete way up to a maximal 
size (nx ~ 100) and by using a continuous variable with an increasing step 
beyond [14]. 



dt 
dCi 



dt 



2.2 Condensation rate 



When the solute long-range diffusion controls the precipitation kinetics, the 
condensation rate is obtained by solving the diffusion problem in the solid 
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solution around a spherical precipitate (cf. appendix B). For a cluster of radius 
r„ x , this condensation rate takes the form [15,16] 

/3n x =47rr nx ^C 1 , (4) 

where D x is the diffusion coefficient of the solute at infinite dilution of the solid 
solution and Q is the mean atomic volume corresponding to one lattice site. We 
checked by studying growth of planar precipitates that precipitation kinetics 
can truly be assumed to be controlled by solute long range diffusion and that 
the diffusion coefficients deduced from kinetic Monte Carlo simulations are in 
good agreement with the experimental ones entering equation 4 [17]. 



2.3 Evaporation rate 



By contrast with the condensation rate, the evaporation rate a nx cannot be 
obtained directly. It has to be deduced from f3 nx using the equilibrium cluster 
size distribution (Eq. 1). 

The evaporation rate is obtained assuming that it is an intrinsic property of 
the cluster and therefore does not depend on the solid solution surrounding 
the cluster 1 . This means that the cluster has enough time to explore all its 
configurations between the arrival and the departure of a solute atom. Thus 
a„ x should not depend on the nominal concentration of the solid solution and 
can be obtained by considering any undersaturated solid solution of nominal 
concentration x x . Such a solid solution is stable. Then there should be no 
energy dissipation. This involves that all fluxes J nx ^ nx +i equal zero. Using 
equation 3, one obtains 



— Aixl^x) 



0, Cnx(4) (5) 



Cnx+lf^x) 



where overlined quantities are evaluated in the solid solution at equilibrium. 
In particular, the equilibrium cluster size distribution is given by equation 1. 
This finally leads to the following expression of the evaporation rate 

Dx 

a nx+ i = Anr nx — exp [(G„ x+ i - G„ x - Gi)/kT\. (6) 



1 Sometimes, one assumes instead that a constraint equilibrium exists for the clus- 
ters in the supersaturated solid solution and that the expression 1 holds even for 
concentrations higher than the solubility limit . Katz and Wiedersich [18] pointed 
that this assumption can lead to the same expression of the evaporation rate as the 
one we used. In particular, this is true when the condensation rate varies linearly 
with the monomer concentration which is the case here. 
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As the condensation rate varies linearly with the monomer concentration, the 
contribution of the effective chemical potential cancels out in the expression 5 
of a nx : we recover our starting assumption. Indeed, all informations concerning 
the solid solution are contained in this chemical potential. Therefore a nx only 
depends on the clusters free energies and not on the overall state of the gas 
of clusters. Especially, the evaporation rate is independent of the nucleation 
free energy, i.e. the free energy gained when nucleating a precipitate out of 
the solid solution: such concepts are useful in the classical nucleation theory 
but play no role in cluster dynamics. This can be seen by dividing the cluster 
formation free energy into a volume and an interface contributions [7], 

G nx - 2n x /i(4 ) = 4n x AG nuc (x° x ) + (36nn x 2 ) 1/3 a 2 a nx , (7) 

AG nuc (x x ) being the nucleation free energy of a solid solution of nominal 
concentration x x , o~ nx the interface free energy of a cluster containing nx solute 
atoms and a the lattice parameter. This leads to the following expression of 
the evaporation rate 



Dx 

n (8) 



exp [(ZQ^a 2 ((n x + 1) 2/ V„ X+1 - n x 2/ V nx - a ± ) /kT] . 



Looking at the expression 4 of the condensation rate and the expression 8 
of the evaporation rate, one sees that the only parameters needed by cluster 
dynamics are the diffusion coefficient and the interface free energy. There is no 
need to know the nucleation free energy 2 in opposition to other mesoscopic 
models based on classical nucleation theory following Langer and Schwartz 
approach [19] as modified by Kampmann and Wagner [1,20-22]. In cluster 
dynamics, thermodynamics of the solid solution is described thanks to a lattice 
gas model and therefore the nucleation free energy results from this description 
and is not an input of the modeling. 



2.4 Definition of precipitates 



The master equation which describes cluster dynamics yields the full cluster 
distribution as a function of time. In order to follow the precipitation kinetics 
and make the link with experimental results, it is convenient to define mean 
values such as the precipitates mean size or their density so as to characterize 

2 Some cluster dynamics studies [12] makes the evaporation rate depend on the nu- 
cleation free energy because a constrained equilibrium of the cluster size distribution 
has been considered so as to reproduce the instantaneous monomer concentration 
C\ and not the total solid solution concentration x x as it should be. 
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kinetics by the evolution of these quantities. Therefore, the solid solution has 
to be discriminated from the precipitates: one needs to define an arbitrary 
threshold size n x above which clusters represent stable precipitates and below 
which they represent fluctuations in the solid solution. One natural choice for 
this threshold size is the initial critical size given by the classical nucleation 
theory, 



This requires to calculate for a given temperature and a given composition 
of the solid solution the nucleation free energy. This cannot be directly done 
with the lattice gas model used by cluster dynamics as the effective chemical 
potential /zx(#x) appearing in equation 7 and needed to obtain AG nuc (x x ) is 
implicit. This potential can be deduced from the simulation when a stationary 
cluster size distribution is observed 3 but it is not known a priori. Nevertheless, 
the nucleation free energy can be estimated by other means and in particular 
thanks to different mean field approximations. In reference 7 it was shown that 
the cluster variation method (CVM) [23, 24] in the tetrahedron-octahedron 
approximation leads to a good prediction of AG" luc (a; x ). Therefore the value 
of the threshold size is chosen equal to that of the critical size in the classical 
nucleation theory when CVM is used to calculate the nucleation free energy. 
Thermodynamics, and then the nucleation free energy, deduced from CVM can 
differ from the ones really corresponding to cluster dynamics. Nevertheless, 
this nucleation free energy is only used to calculate a threshold size and as the 
aim of this paper is to compare precipitation kinetics obtained with Monte 
Carlo and cluster dynamics simulations, this is not a real issue. The key point 
is to use the same threshold size when comparing the two modeling techniques. 



3 Cluster dynamics simulations 

In references 7 and 8 one of us developed an atomic kinetic model for Al-Zr 
and Al-Sc systems. Parameters were deduced from experimental data (Zr and 
Sc solubility limits in Al and diffusion coefficients) as well as from ab-initio 
calculations (Ah>Zr and ALjSc free energy of formation). So as to compare 
the precipitation kinetics obtained from cluster dynamics with kinetic Monte 
Carlo simulations we have to deduce from this atomic model the interface free 
energies between the precipitates and the aluminum matrix. As for diffusion 
coefficients, we use in the cluster dynamics simulations the experimental values 
known for Zr and Sc impurity diffusion in aluminum [25-27]: 



3 The effective chemical potential is then linked to the instantaneous monomer 
concentration by the relation \i = \G\ + fcTln(Ci)] /2. 
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£> Zr * = 728 x 10~ 4 exp(-2.51 eV/kT) m 2 .s~\ 
D Sc * = 5.31 x 10~ 4 exp(-1.79 eV/kT) m 2 .s _1 . 

As shown in Ref. 7,8, the diffusion model used in kinetic Monte Carlo does 
reproduce the above values. 

3. 1 Interface free energy 

The interface free energy can be deduced from the atomic model by computing 
the free energies corresponding to planar interfaces in the three most dense 
packing orientations {111}, {110} and {100} within the Bragg- Williams ap- 
proximation and then by using the Wulff construction to obtain an isotropic 
average interface free energy a [7]. Such an interface free energy depends on 
the temperature for two reasons: atomic interactions used to obtain a vary 
with temperature and there is a configurational entropy contribution due to 
the relaxation of the interface. On the other hand it does not depend on the 
size of the cluster and corresponds to the asymptotic limit of cr nx . A direct 
calculation considering the partition functions of small clusters shows that <x nx 
slightly deviates from this asymptotic value at small sizes [7]. In the following, 
we test effects on cluster dynamics results due to this slight dependence of 
the interface free energy with the cluster size. Some authors [9,11,28] already 
used an interface free energy depending on the cluster size so as to model 
copper precipitation in iron by cluster dynamics. In these studies, the size 
dependence was introduced to take into account the change with size of the 
precipitate structure. This is not the purpose here as precipitates keep their 
Ll 2 structure whatever their size. We want to show how sensitive are cluster 
dynamics simulations to the interface free energy and why one needs to go 
beyond a model using a constant parameter. 

In this purpose, for clusters containing no more than 9 solute atoms, the in- 
terface free energy is computed using Eq. 20 of Ref. 7. For clusters of size 
nx > 10, the interface free energy is obtained using an extension of the capil- 
lary approximation, 

o- nx = a (l + c n x 1/3 + d n x 2/3 j , (10) 

where c and d respectively correspond to the line and point contributions 
[29]. The asymptotic value a used is the one previously calculated thanks 
to the Bragg- Williams approximation and the Wulff construction [7] whereas 
coefficients c and d are obtained by a least square fit of the expression 20 of 
Ref. 7 for sizes 5 > nx > 9. The curvature correction c obtained with such a 
procedure is found to be negative in agreement with classical models [30,31]. 
Instead of using such an approximate procedure to obtain the interface free 
energy, this can be done more precisely by sampling thermodynamic averages 
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with a Monte Carlo algorithm so as to compute the free energy difference 
between a cluster of size n and one of size n + 1 at a given temperature [29]. 
Interface free energies obtained by both methods are in good agreement [32]. 

As seen in figure 1, the difference between this interface free energy <7„ x de- 
pending on the cluster size and its asymptotic limit a is small. For instance, for 
the interface between Al 3 Sc and the solid solution at T = 500°C, the maximal 
difference is obtained for sizes % ~ 10 and is less than 5%. 

We run some cluster dynamics simulations with the interface free energy de- 
pending on the cluster size as given by Eq. 20 of Ref. 7 and Eq. 10 and 
some other simulations with a constant interface free energy taken equal to 
the asymptotic limit a. The comparison with kinetic Monte Carlo simulations 
(Fig. 2) shows that one obtains a better agreement when taking into account 
the size dependence: the number of precipitates as well as their mean size 
predicted by both simulation methods are then quite close. Although the size 
dependence of the interface free energy is small, it definitely improves the 
ability of cluster dynamics to well reproduce kinetic Monte Carlo results. This 
is true mainly for the precipitate density, the effect on precipitate mean size 
being less pronounced, especially when the supersaturation is low. The density 
is more sensitive than the mean size to the distribution of small clusters. As 
these clusters are the most affected by the variation with size of the interface 
energy, this explains why the greatest effect observed is for the precipitate 
density. 




Fig. 1. Variation with the cluster size ns c of the interface free energy between the 
solid solution and AI3SC at T = 500°C. Symbols correspond to a ngc as given by the 
direct calculation of the cluster formation free energy (Eq. 20 of Ref. 7) and lines 
to the capillary approximation (Eq. 10) as well as to its asymptotic limit. 
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I (s) I (s) 

(a) x° Sc = 0.4 at.% (n* Sc = 16) (b) x° Sc = 0.75 at.% (n* Sc = 9) 

Fig. 2. Precipitation kinetics of two supersaturated solid solutions of nominal con- 
centration x^ c = 0.4 at.% and x^ c = 0.75 at.% at T = 500°C: evolution with time 
of the number N p of precipitates in the simulation box (normalized by the number 
of lattice sites N s ) and of their average size (rasc)p- Full lines correspond to cluster 
dynamics simulations with an interface free energy depending on the cluster size 
(Eq. 20 of Ref. 7 and 10), dotted lines to simulations with a constant interface free 
energy and symbols to kinetic Monte Carlo simulations 

In the following, all cluster dynamics simulations will use this size dependent 
interface free energy as input parameter. 



3.2 Precipitation kinetics 

Comparing, at a given temperature, the precipitation kinetics obtained with 
Monte Carlo and cluster dynamics simulations for different nominal concen- 
trations of the solid solution (Fig. 3), one sees that cluster dynamics manages 
to reproduce the variations of the precipitate density as well as the variations 
of their mean size. For low supersaturations, the agreement is really good 
whereas for higher supersaturations there is a delay: cluster dynamics is too 
slow compared to kinetic Monte Carlo simulations by a factor of roughly 2. 
Nevertheless, in all cases, the prediction of the precipitate maximal density 
at the transition between the growth and coarsening stages is correct. Cluster 
dynamics technique manages to reproduce the variations over three orders of 
magnitude of this maximal density in the range of concentrations considered 
in Al-Sc solutions at T = 450°C (Fig. 3). 
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Fixing now the nominal concentration of the solid solution and varying the 
temperature, the comparison of results obtained with both simulation tech- 
niques leads to the same conclusions (Fig. 4). Cluster dynamics results still 
appear to be too slow at high supersaturations, but except this time delay, 
variations of precipitate density and of their mean size are well reproduced. 





^ — 3 — ° °- 1 ' 1 1 1 

0.0001 0.001 0.01 0.1 1 

f(S) 



Fig. 3. Precipitation kinetics for different nominal concentrations of a supersatu- 
rated aluminum solid solution at T = 450°C. Lines correspond to cluster dynamics 
simulations and symbols to kinetic Monte Carlo simulations. 
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0.01 0.1 1 10 100 1000 0.001 0.01 0.1 1 10 100 1000 



«(s) t(s) 

(a) x\ r = 1 at.% (b) x\ x = 1.25 at.% 

Fig. 4. Precipitation kinetics for different temperatures of a supersaturated solid 

solution of nominal concentration x? r = 1. at.% and x? r = 1.25 at.% 
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3.3 Cluster size distribution 



Cluster dynamics lead to a stationary cluster size distribution in the solid 
solution, i.e. for clusters smaller than a critical size (Fig. 5), as predicted 
by classical nucleation theory. This stationary distribution corresponds to a 
constraint equilibrium existing for sub-critical clusters but nothing has been 
supposed concerning clusters bigger than the critical one. Using the capillary 
approximation to estimate the cluster formation free energy (Eq. 7), classical 
nucleation theory assumes this stationary distribution to be given by equa- 
tion 1. Obviously, in this expression, the interface free energy has to be taken 
equal to the one used for cluster dynamics simulations and therefore we con- 
sider the size dependent interface free energy given by Eq. 20 of Ref. 7 and 
Eq.10. As for the nucleation free energy AG nuc (x^), we already pointed that 
several thermodynamic approximations can be used to estimate it. One usu- 
ally considers ideal or regular solid solution models to calculate AG nuc (xx) but 
it was shown in reference 7 that at least CVM in the tetrahedron-octahedron 
approximation has to be used in order to obtain a good agreement between 
classical nucleation theory predictions and kinetic Monte Carlo simulations. 
The comparison of the stationary cluster size distributions in the solid solu- 
tion obtained from cluster dynamics simulations with the ones predicted by 
classical nucleation theory leads to the same conclusion: one has to use CVM 
to calculate the nucleation free energy in order to get the right distribution 




10 20 30 40 



Fig. 5. Cluster size distribution of an aluminum solid solution of nominal concentra- 
tion Xg c = 0.2 at.% at T = 450°C. Symbols correspond to cluster dynamics simula- 
tions and lines to prediction of classical nucleation theory with different mean-field 
approximations of the nucleation free energy. 
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(Fig. 5) 4 . This result can appear quite surprising: although no information 
concerning the nucleation free energy has been used to run cluster dynamics 
simulations, this method leads to the stationary size distribution correspond- 
ing to the CVM. It seems that the lattice gas thermodynamic description 
used by the cluster dynamics manages to take into account short range order 
effects which are considered within CVM, in opposition to ideal and regular 
solid solution approximations. 



4 Effect of the overlap of diffusion fields 

As seen in the previous section, discrepancies between cluster dynamics and ki- 
netic Monte Carlo simulations appear for high supersaturations, around 1 at.% 
at T = 450°C for instance. Although the maximal precipitate density and its 
variations with temperature and supersaturation are correctly predicted, the 
precipitation kinetics as given by cluster dynamics are too slow. One possible 
origin of such a discrepancy might be in the expression for the condensation 
rate /3„ x : the later is computed solving, under stationary conditions, the diffu- 
sion equation around the cluster isolated in an infinite medium. This assump- 
tion looks reasonable for dilute solutions but its validity is not so clear for high 
supersaturations when clusters are quite close from each other: overlap of dif- 
fusion fields might alter the rate of solute impingement at clusters. Therefore, 
in order to see if a better agreement can be obtained between cluster dynam- 
ics and kinetic Monte Carlo simulations, it is worth trying to describe kinetic 
effects arising from interactions between the diffusion fields created around 
each precipitate and to incorporate those in cluster dynamics simulations. 

In this purpose, we follow the approach proposed by Brailsford et al. [34] and 
already used by Smetniansky-De-Grande and Barbu in their modeling of the 
precipitation kinetics in Fe-Cu alloys [35,36]. When calculating the monomer 
condensation rate on a given cluster, we assume that at distances greater 
than r ext , one half of the mean distance between precipitates, monomers feel 
influence of the cluster through its diffusion field as well as the one arising 
from other clusters through an effective medium characterized by a parameter 
k. The rate at which this effective medium absorbs monomers is given by 
D x k 2 (Ci — Ci q ), C\ and C{ q being respectively the instantaneous and the 
equilibrium monomer concentrations. This leads to the following expression 
of the condensation rate, depending on whether the radius r„ x of the cluster 



4 The difference between the distribution given by equation 1 and the observed 
stationary distribution during cluster dynamics simulations for size equal or slightly 
smaller than the critical one can be attributed to the Zeldovitch factor. Indeed, for 
the critical size, there should be a factor 1/2 between the two distributions [33]. 
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is smaller or greater than r ext (cf. appendix B), 

^ = 4 "^ C ' l + ^ r -r„ x) - r "* Sr " (lla) 



A, x = 47T^r nx C 1 (l + kr„) , r„ x > r"'. (1110 



Comparing this expression with the previous one (Eq. 4) where precipitates 
were assumed to be isolated, we see that the consideration of interactions 
between the different cluster diffusion fields through the introduction of this 
effective medium leads to a greater condensation rate. 

So as to ensure that cluster dynamics simulations converge to equilibrium, 
i.e. a two-phases system consisting of AI3X precipitates embedded in a solid 
solution of composition x^, evaporation rates need to be multiplied by the 
same factor as the condensation rates. They are now given by 



1 + kr ext 
1 + k(r ext - r„ x _i y 



««x = o£ X i I UZ^t, ~ V' r «x-i < r ex \ (12a) 



« nx = o£ x (l + fcr nx _i) , r nx _i > r ext , (12b) 

where a° x is the evaporation rate of an isolated cluster as given by equation 
8. Through the parameters k and r ext , the evaporation rates now depend on 
the solid solution considered. 



As the effective medium we introduced has to be equivalent to the cluster 
assembly, the parameter k is calculated so as to ensure that the quantity of 
monomers absorbed by this effective medium is equal to the quantity con- 
sumed by all clusters, 

D x k 2 (d - Cl q ) = 2J^ 2 + ]T J n „ +1 . (13) 

n x >2 

This equation has to be solved self-consistently at each step of cluster dynam- 
ics simulations. Starting from the value of k corresponding to the previous 
iteration or from for the first step, equations 11 and 12 allow to calculate 
the cluster condensation and evaporation rates and thus the different fluxes. 
Using equation 13 a new value of the parameter k can be computed. To do 
so, the equilibrium monomer concentration is required. It is simply given 
by the equilibrium cluster size distribution 1 when the nucleation free energy 
is null, 

C 1 e9 = exp(-(367r) 1/ V 1 /A;T). (14) 

Precipitation kinetics obtained with cluster dynamics when considering in- 
teractions between the cluster diffusion fields are not really different from 
the ones previously obtained without such an interaction. At low supersatu- 
rations, no difference is observed whereas at higher supersaturations (Fig. 6) 
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this interaction slightly speeds up the precipitation kinetics. Nevertheless, this 
does not really improve the agreement with kinetic Monte Carlo simulations 
and kinetics obtained with cluster dynamics still appear too slow. Therefore 
the discrepancies appearing at high supersaturations do not arise from the 
assumption of isolated clusters used to calculate the condensation rate but 
should have another origin. We checked too that the use of stationary condi- 
tions to compute the condensation rate is correct by using the whole solution, 
i.e. incorporating the transition state, of the diffusion equations [15] but this 
does not lead to a better agreement. Another possible origin of the discrepancy 
could be the use of a single parameter, their size, to describe clusters which 
involves considering clusters as spherical particles. This assumption might be 
too crude when computing the condensation rate but we did not manage to 
check this point as this spherical approximation is needed to solve the long- 
range diffusion equations. 

As the correction taking into account the overlap of cluster diffusion fields 
does not really change cluster dynamics simulations, we will not consider it 
in the following section when comparing simulated kinetics with experimental 




• 1 1 1 1 

1 10 100 1000 10000 

t(s) 

Fig. 6. Precipitation kinetics of a supersaturated solid solution of nominal concentra- 
tion x^y. = 1 at.% at T = 450°C. Lines correspond to cluster dynamics simulations 
with or without interaction between the diffusion fields created around each cluster 
and symbols to kinetic Monte Carlo simulations. 
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data: the condensation and evaporation rates will be those of isolated clusters 
as given by equations 4 and 8. 



5 Comparison with experimental data 

Unlike kinetic Monte Carlo which can only simulate rapid processes (early 
stages of precipitation, high supersaturations), cluster dynamics can be used 
for low supersaturations and long annealing times. This allows a direct compar- 
ison with experimental studies. The precipitation kinetics in Al-Zr [1-3,37-40] 
as well as in Al-Sc system [4-6,22,41-46] has been studied by several groups. 
However, for the Al-Zr system, experimental data can hardly be directly com- 
pared to cluster dynamics simulations. Indeed, one observes two kinds of pre- 
cipitates, some spherical ones and some rod-like shaped ones. Moreover, for 
low supersaturations, nucleation occurs heterogeneously and, Al-Zr being peri- 
tectic, high supersaturations are unreachable. Another difficulty arises from 
the fact that the Ll 2 structure is unstable and that the stable DO23 struc- 
ture of Al 3 Zr precipitates appears for long enough annealing times. On the 
contrary, Al-Sc system appears to be a good candidate to compare simulated 
precipitation kinetics with experimentally observed ones as experimental data 
allowing to quantify the homogeneous precipitation of the Ll 2 structure of 
Al 3 Sc exists. 

Novotny and Ardell [6] as well as Marquis et al. [5,44,45] observed the coarsen- 
ing behavior of different Al-Sc alloys. For the higher supersaturations studied 
(a;g c = 0.18 at.%), the authors concluded that precipitation occurs homo- 
geneously, at least for temperatures not too high (T < 350°C). Therefore, 
precipitation kinetics of the same alloy can be simulated using cluster dynam- 
ics. For the three different temperatures studied (T = 300, 350, and 400°C), 
it appears that the cluster dynamics equation manage to reproduce the varia- 
tion with time of the mean precipitate radius (Fig. 7): the agreement is really 
good with Novotny's data and reasonable with Marquis' ones. Notice that the 
parameters values used in the cluster dynamics were deduced from the atomic 
diffusion model used in the kinetic Monte Carlo simulations of the short an- 
nealing times of the same Al-Sc alloy maintained at the same temperatures. 
Cluster dynamics, with a single set of parameters, thus reproduces atomic 
simulations at short times and gives a safe extrapolation thereof to the range 
of annealing times that can be compared with experimental data. It should 
be noticed too that cluster dynamics manages to catch the variation with the 
annealing temperature of the coarsening kinetics: the experimental observed 
speed up of the cluster growth rate by a factor ~ 10 when going from T = 300 
to 400°C is well predicted. The good agreement found for the lower tempera- 
tures (T = 300 or 350°C) depreciates a little bit at higher temperatures. This 
is due to the fact that, experimentally, for a given nominal concentration one 
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Fig. 7. Mean precipitate radius as a function of the aging time for a solid solution 
of composition zg c = 0.18 at.% at different temperature (T = 300, 350, and 400°C) 
as given by cluster dynamics and compared to experimental data [5,6,44,45] and to 
kinetic Monte Carlo results [7,8]. The cutoff radius used for cluster dynamics and 
KMC is r£ ~ 0.75nm (n£ = 27). 
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reduces the supersaturation by increasing the temperature and thus favors 
heterogeneous compared to homogeneous precipitation. There are then less 
precipitates, these ones being larger. 

In figure 8, we compare the precipitate size distribution obtained from cluster 
dynamics with the experimentally ones measured by Novotny and Ardell [6]. 
For the different annealing times studied the two normalized distributions are 
close. Actually, cluster dynamics leads to a normalized size distribution not 






Fig. 8. Normalized precipitate size distributions g(r) for an Al-0.18 at.%Sc solid 
solution at T = 350° C obtained with cluster dynamics simulations and compared 
to experimental data [6] as well as to LSW distribution. 
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really different from its asymptotic limit, the Lifshitz-Slyozov- Wagner (LSW) 
distribution. Novotny and Ardell already concluded to a good agreement be- 
tween experimental size distributions and the LSW distribution, which cor- 
responds to the same agreement found here with the distributions obtained 
from cluster dynamics. 

Hyland [4], using TEM, measured the variation with time of the precipi- 
tate density for two different annealing temperatures in an alloy containing 
0.11 at.% Sc. It was not possible to reproduce satisfactorily these experimen- 
tal data with cluster dynamics. Indeed the simulated densities are higher by 
at least one order of magnitude than the experimental ones. The number of 
precipitates experimentally observed depends on the spatial resolution of the 
observation technique. This can be reproduced by imposing in the simulations 
a cutoff radius below which the clusters are ignored, even though they might 
be supercritical. The densities calculated from the simulations are really sen- 
sitive to the cutoff radius used to count precipitates. It is possible to obtain 
lower densities by increasing this cutoff radius but, then, the time scale of the 
simulations completely differs from the experimental one. The origin of this 
discrepancy remains an open question. 



6 Conclusions 

Precipitation kinetics of Al 3 Zr and Al 3 Sc in aluminum supersaturated solid 
solution has been modeled using the cluster dynamics technique. The only 
input parameters are the solute diffusion coefficients and the precipitate / solid 
solution interface free energies. In contrast with other mesoscopic modeling 
techniques based on Kampmann and Wagner approach, cluster dynamics in 
its strict sense does not require a definition of the nucleation free energy: 
this is not an input parameter of the model but rather a by-product of the 
cluster gas thermodynamic underlying the cluster dynamics equations. We 
showed that this description leads to a stationary cluster size distribution 
in the metastable solid solution corresponding to the one predicted by the 
classical nucleation theory provided CVM is used to calculate the nucleation 
free energy, in agreement with previous study. 

Cluster dynamics is shown to be very sensitive to the interface free energy: the 
size dependence of the interface free energy, which is revealed by the direct 
computation of the free energy of small clusters, is crucial in making the 
cluster dynamics calculations reproduce kinetic Monte Carlo simulations based 
on the same atomistic parameters. When doing so, the precipitation kinetics 
obtained agree well with kinetic Monte Carlo simulations: the variations of the 
precipitate density and of their mean size are well reproduced. Nevertheless, 
for high supersaturations, kinetics modeled by cluster dynamics appear to 
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be slightly too slow. We showed that this discrepancy cannot be solved by an 
improved kinetic model taking into account the overlap of the cluster diffusion 
fields. 

The cluster dynamics technique as described manages to well reproduce ex- 
perimental data, despite the fact that the time scales in real experiments 
are several order of magnitude larger than those in kinetic Monte Carlo. The 
present study therefore validates, in a quantitative manner, the proposed mul- 
tiscale scheme for homogeneous precipitation kinetics: atomic scale parameters 
are obtained from available experimental quantities and ab-initio calculations; 
mesoscale parameters are deduced from the latter. Cluster dynamics as used in 
this work, i.e. in its strict sense, then appears as the good tool to extrapolate 
kinetic Monte Carlo simulations, which reveal full details of the nucleation 
process, to a wider range of annealing times, of temperatures and of supersat- 
urations. It provides an effective method to accomplish quantitative transition 
of scales, at least in the low solubility systems which have been considered in 
this work. 



A Equilibrium cluster size distribution 

In cluster dynamics, the size distribution of the clusters in a solid solution at 
equilibrium is used to deduce the monomer evaporation rates from the con- 
densation rates. This distribution is needed by classical nucleation theory too 
so as to estimate the cluster size distribution in the metastable solid solution 
and to deduce from it the steady-state nucleation rate. In the following, we 
show how this distribution can be obtained for an undersaturated solid solu- 
tion, as needed by cluster dynamics, the classical nucleation theory extending 
the validity of this distribution to supersaturated solid solutions. 

The cluster size distribution in a solid solution at equilibrium, can be estab- 
lished quite easily as long as the solid solution is dilute [47,48]. Considering 
an assembly composed of N nx clusters containing nx solute atoms and of free 
energy G nx , the total free energy of the system can be written 

oo 

N S G = N n x G nx + kTlog(W), (A.l) 

n x =l 

where N s is the number of sites that can be occupied by a cluster and W is the 
number of different configurations accessible to the cluster assembly. Assuming 
that each cluster whatever its size lies only on one site and neglecting around 
each cluster all excluded sites which cannot be occupied by any other cluster, 
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this number is simply given by 



N ' 

W = 7 — -r . (A.2) 

y ly s Z^n x =l "x^- lln x =l " 



"n x - 



The concentration of each size class being C„ x = N nx /N s , Stirling formula 
leads to the following estimation for the total free energy, 

oo oo 

G= E C nx G nx + kT ]T C„ x log(C nx ) 

(A.3) 

+ A;T ( 1 - 2 Cn x I log ( 1 - £ C nx ' 

\ "x=l / \ "x=l 

The equilibrium concentrations (7 nx are obtained by minimizing this free en- 
ergy while imposing that the solute concentration equals a fixed value x x : 

oo 

Y: n x C nx = 4. (A.4) 

n x =l 

This can be done quite easily by considering the grand canonic free energy 

oo 

A=G+ (1-2 £ n x C nx V, (A.5) 

n x =l 

where we have introduced a Lagrange multiplier /i to ensure the constraint 
A.4. /i = (fix — / i A)/2 is nothing else than the effective chemical potential, 
i.e. half the difference between the chemical potentials of the solvent A and 
the solute B. The minimization of this grand canonic free energy leads to the 
following solution for the cluster size distribution 

_^ = e ^(-^*i). (A,, 

AG„ X = G nx — 2nxH is the cluster formation free energy relative to the solid 
solution. It can be estimated using the capillary approximation. As we made 
the assumption of a dilute solid solution to estimate the total free energy G, 
the sum appearing in the left hand side of equation A.6 can be neglected, so 
the equilibrium cluster size distribution is simply given by 

C , nx =exp(-AG nx /A;T). (A.7) 



B Condensation rate 



Monomer condensation rates on clusters are the main input of cluster dynam- 
ics. Assuming that the precipitation kinetics is controlled by the solute long 
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range diffusion, the condensation rates are obtained by solving the diffusion 
equation for the solute in the solvent matrix to a spherical cluster with the 
following boundary conditions: the solute concentration is zero at the cluster 
(in order to describe solute trapping at the cluster) and far from this cluster 
the monomer concentration is equal to that corresponding to cluster dynamics 
simulation. Two cases can be distinguished, depending on whether the cluster 
is assumed to be isolated or embedded in an effective medium which simulates 
the trapping to other clusters, i.e. the effect of the overlap of diffusion fields. 



B. 1 Isolated precipitate 



The rate at which monomers condense on an isolated cluster of radius r„ x is 
obtained by solving the stationary diffusion equation 

D x AC 1 (r)=0, (B.l) 

with the following boundary conditions: Ci(r nx ) = at the cluster- matrix 
interface 5 and C\(r — > oo) = far away from the cluster. The monomer 
concentration profile around the precipitate is thus 

d(r) = C 1 °° r ~ rrax . (B.2) 
r 

The condensation rate is given by the integral over the cluster surface of the 
incoming flux of monomers, 



3 =-A,r 



'"X 



where f2 is the mean atomic volume corresponding to one lattice site. 



(B.3) 



B.2 Precipitate in an absorbing medium 



In order to account for the overlap of the diffusion fields around the many 
clusters in the matrix, one can compute the diffusion flux to one specific cluster 
embedded in an effective absorbing medium. The latter extends from infinity 

5 Actually, at the interface between the cluster and the matrix, the monomer con- 
centration is equal to the equilibrium one corrected from Gibbs-Thomson effect. 
Taking this concentration equal to is equivalent to assume that the cluster does 
not emit any monomer and allows one to directly obtain the condensation rate. 
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to a distance r ext of the cluster, equal to the mean inter-cluster distance. This 
effective medium is characterized by a constant k giving the rate at which 
monomers are absorbed, D x k 2 (Ci(r) — C{ q ). Therefore, for precipitates of 
radius smaller than r ext , the condensation rate is obtained by solving the set 
of equations 

D x Ad(r) = 0, r < r ext (B.4a) 
D x ACi(r) - D x k 2 (d(r) - C{ q ) = 0, r> r ext . (B.4b) 

The same boundary conditions apply as previously on the monomer concentra- 
tion with the additional condition that the diffusion flux has to be continuous 
at r = r ext . The resolution of this diffusion problem gives the following con- 
densation rate, 

D x noo 1 + kr ext 
^ = 47r 7T^ l + k(r^-r nx Y (R5) 

For precipitates of size r n > r ext one has simply to solve the equation 

£> x ACi(r) - D x k 2 (d(r) - C{ q ) = 0, (B.6) 

still with the boundary conditions Ci(r„ x ) = and C\{r — > oo) = . One 
obtains for the condensation rate 

Pn x = 4vr^r rix C 1 00 (l + kr nx ). (B.7) 
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